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ABSTRACT 

We discuss and compare two alternative models for the two-point angular correlation 
function of galaxies detected through the sub-millimetre emission using the Herschel 
Space Observatory. The first, now-standard Halo Model, which represents the angular 
correlations as arising from one-halo and two-halo contributions, is flexible but com- 
plex and rather unwieldy. The second model is based on a much simpler approach: we 
incorporate a fitting function method to estimate the matter correlation function with 
approximate model of the bias inferred from the estimated redshift distribution to find 
the galaxy angular correlation function. We find that both models give a good account 
of the shape of the correlation functions obtained from published preliminary studies 
of the HerMES and H- ATLAS surveys performed using Herschel, and yield consistent 
estimates of the minimum halo mass within which the sub-millimetre galaxies must 
reside. We note also that both models predict an inflection in the correlation function 
at intermediate angular scales, so the presence of the feature in the measured correla- 
tion function does not unambiguously indicate the presence of intra-halo correlations. 
The primary barrier to more detailed interpretation of these clustering measurements 
lies in the substantial uncertainty surrounding the redshift distribution of the sources. 



1 INTRODUCTION 



. Current theoretical models predict that galaxies form and 
' evolve in cold dark matter (CDM) halos. Galaxies conse- 
, quently tend to trace the distribution of mass, although the 
manner in which they do this may be biased (Kaisei| [l984l : 
LBardccn et al. 1986; Coles 1993; Mo & White 1996). Inprin- 
, ciple, therefore, once the bias is allowed for, it is possible to 
' use measurements of the clustering of galaxies to determine 
the clustering properties of the dark matter, especially if 
measurements can be made as a function of redshift. Because 
clustering evolution is sensitive to the parameters underly- 
ing the background cosmological model, such observations 
can provide another (independent) test of the concordance 
cosmological model; see, e.g. IColed (|2005l ). In addition, clus- 
tering observations can be used to constrain properties of 
the galaxies themselves, such as the minimum halo mass 
within which they reside, which may yield clues about the 
processes of galaxy formation and evolution. 

A steadily increasing number of surveys of large-scale 
galaxy clustering are now becoming available. In the opti- 
cal wavebands there are projects s uch as the UKIDSS Ultra 
Deep S urvev (iHartlev et al.ll2010l) and the SPSS R edshift 
Survey (jConnollv et al.ll20ld : IRoss fc Brunneij 12009) which 
are being used to extract information on clustering as a 
function of redshift. The Herschel Space Observatory was 
launched in 2009 and is the only space observatory to cover 
a spectral range from the far infrared to sub-millimetre and 
therefore provides a new and unique window through which 



to study high-redshift galaxy clustering. Two surveys of par - 
ticular interest to this article, H erMES jOliver et al.ll2O10l ) 
and H- ATLAS l|Eales et aL 2010l). have already released an- 
gular correlation results ( Coorav et al.|[2010l : iMaddox et al.l 
l20ld ) from their Science Demonstration Phase (SDP) which 
we will discuss further later, so it is timely to raise the issue 
of modelling sub-millimetre galaxy clustering. 

The so-called Halo Model has been used extensively 
over the past few years in modelling galaxy clustering in 
a variety of contexts. The Halo Model combines approxima- 
tions of the dark matter profile within individual halos, the 
mass function, and bias models to estimate the correlation 
function for given cosmological parameters and characteris- 
tic halo masses. It has been shown to provide accurate and 
reliable predictions of clustering measurements, at the price 
of a certain degree of complexity and modelling freedom. 

The main focus of this paper is a comparison of the Halo 
Model and a fitting function method which was initially in- 
troduce d bvlHam ilton ct al. (1991), and subsequently dcvel - 
oped bv lPeacock fc Doddd (119941) . I Jain. Mo fc Whit j (| 19951 ) 
and especially Peacock fc Doddj 1 19961) . to estimate the 
matter corr elation function. Work b y Matarrese et al.l (Il997l ) 
followed bv lMoscardini et al.1 (|l998l ) and lColes et all (| 19981 ) . 
showed how to incorporate this idea into a technique for 
providing detailed predictions of angular correlation in high- 
redshift galaxy surveys. In this paper we compare the pre- 
dictions of this much simpler approach with results from the 
Halo Model. 
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The paper is organized as follows. In the next section, 
we outline the methodology for the fitting function method. 
Following that in Section O the models are compared to 
recent Herschel data and the best-fitting values of the free 
parameters are extracted and compared. In Section |4j we 
conclude with a summary. 



2 MODELS OF GALAXY CLUSTERING 

A model for the angular correlation function of galaxies basi- 
cally involves a model for the three-dimensional clustering of 
galaxies as a function of epoch (i.e. redshift z) combined with 
a machinery (derived from a model of the background cos- 
mology) for projecting this information down the observer's 
light cone. 



2.1 Angular correlation function 

The observed angular correlation function, cjobs iP) , is calcu- 
lated from the spatial power spectrum, A^{k,z), using: 



uJobs{9)= G{z)Af\z) / ^{r{u,e,z),z)dudz 



(1) 



where M{z) is the normalized distribution of galaxies over 
redshift, G{z) — {dx/dz)~^ where x{z) is the comov- 
ing distance to redshift z, the separation r{u,6,z) = 
^Jv? + x^{z)9'^ Mpc^~"^, and ^{r,z) is the two-point 
(galaxy-galaxy) correlation function where, 



ar,z)= / Alik,z) 



sin(fcr) dk 
kr k 



(2) 



The galaxy power spectrum can be approximated using dif- 
ferent methods such as those described below. 



2.2 Power Spectrum: Halo Model 

The Halo Model estimate for the galaxy power spectrum 
represents it as arising from two distinct components: 



Al{k,z) = Al^{k,z) + Al„{k,z), 



(3) 



The first term represents contributions from galaxies resid- 
ing within the same dark matter halo, while the second is 
generated by contributions from galaxies in separate halos. 
The one-halo term models small-scale highly non-linear clus- 
tering using a recipe that involves a number of different com- 
ponents, including the dark matter profile of the halos, the 
mass function of halos, and the halo occupation distribu- 
tion (HOD); the two-halo term depends on the underlying 
matter power-spectrum as well as the bias of galaxy clus- 
tering relative to mass clustering. In this paper the dark 
matter halo density profile is approximated using the NFW 
profile suggested by iNavarroT Frenk fc Whit3 (|l997l ). The 
mass function is effectively the number density of halos at 
a given mass and fo r this we use the approximation found 
bv lsheth fc TormenI (|l999l '). The Halo Occupation Distribu- 
tion (HOD) is a simple way of relating the distribution of 
galaxies from the dark matter distribution, by estimating 
how many galaxies are expected within a dark matter halo 



of given mass. This is split into two estimates - the number 
of central galaxies (which would either be 1 or 0) and the 
number of sat ellite galaxies; a summary of the procedure can 
be found in , Ross fc Brunnej 12009,1. There are extensions to 
the model (e.g. Gil-Marfn. Jimenez, fc Verde 2010l') but we 
just consider the basic version here. lCoorav fc ShethI (|2002D 
provide a detailed review of the Halo Model, so we will not 
go into any further details. 



2.3 Power Spectrum: Fitting Function 

The linear matter power spectrum can be calculated as: 



Alin(fclin 



27r2 



-T\knn)Dl< 



(4) 



where the transfer function, r(fciin), and the growing mode 
of linear p erturbations, D+(z), ca n be approximated as de- 
scribed bv lEisenstein fc Hul (|l999l ). 

The non-linear power spectrum (A,^i) is known for ex- 
treme values of the wavenumber fc (i.e. at small fc, A^i ~ 
Af;jj, and at large fc, A^i ~ Afj^^), the latter derived using 
the so-called stable clustering ansatz. In the light of these 
two asymptotic regimes a fitting function was proposed to 
extrapolate A^j at all other fc directly from the linear coun- 
terpart: 



Anl (fcnl , z) = J^[Aun {hin , z) 



(5) 



where fcni = fcun [1 + A^i (fcni ,z)]^. IPeacock fc Doddd (|l996h 
found a fit for this for low density Universes with cold dark 
matter and a cosmological constant: 



where. 



1 + B/3x + [Ax] 



l + {[Ax]"g^{z)/\Vx^/'^]y 



1//3 



(6) 



A = 0.482(1 4- neff/3)-°■^*^ 
B = 0.226(H-ncff/3)"^■^^^ 
a = 3.310(H-n,ff/3)~°''**, 
= 0.862(1 + n,ff/3)-"■''*^ 
1/ = 11.55(1 -f-neff/3)-''-'*'^ 



and Hcff = 



din 



din fc] 



— ( -3-Aun(fclin' 



This method can be used to approximate the non-linear 
dark matter power sp ectrum (A^m ). Under very general 
conditions, proved by IColed (|l993l '). this can used to ap- 
proximate the non-linear galaxy power spectrum, Ag, via 
a simple linear relationship of the form: 



Ag(fc,2) = KftAjy^{k,z). 
We 



(7) 



use an app r oxima tion for the bias proposed by 
Moscardini et al. f (| 19981 ). which is based on work by 
Mo fc White! (|l996l ) who suggested the following relation be- 
tween the bias and mass density fluctuations: 



b{z, m) 



1 + 



(8) 



where 5c is the critical linear density required for collapse 
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Figure 1. The normalized distribution of galaxies as a function 
of redshift for data sets from Her MES and H- ATLAS. HerMES 
estimates are from ICoorav et al. and are shown in thin 

lines: S260 > 35mjy is a solid line, S350/S250 > 0.85 is the dotted 
line and S350 /S250 < 0.85 is the dashed line. H- ATLAS estimates 
are from lEales et al" I l l2010al ') and are shown in thick lines: S250 > 
33mjy is the dotted line, and S350 > 35mjy is the solid line. 

and o-fi^ is the linear variance in the smoothed density field. 
The effective bias can then be derived as, 

00 

6cff(z) = J b(z,m)n{z,m)dlnm, (9) 

ln("min) 

where n is the nor malized mass function, given by 
ISheth fc Tormenl(| 19991 ). Note for simplicity we assume here 
that Mmin is constant over redshift; given the larger redshift 
ranges a minimum mass varying with redshift may provide 
an improved fit but this approximation is sufficient for the 
comparison of models in this work. 

Apart from the cosmological parameters, which we con- 
strain to match other observations in this analysis, this fit- 
ting function method has just one free parameter (the min- 
imum halo mass Mmin) whereas the Halo model has three: 
the minimum halo mass, the average mass of a halo with one 
satellite galaxy (Msat ) and the slope of the first moment of 
the satellite galaxy HOD (as). 



3 FITS TO THE DATA 

In this section, the models described above are compared 
with the observed angular correlation from recent Herschel 
surveys and are used to put estimates on the minimum halo 
mass associated with these sources. The models we choose 
use current cosmological parameters fr om the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) jKomatsu et al.ll2009l . 
Clm = 0.274, Qa = 0.726, S^i, = 0.046, h = 0.705, n = 0.96, 
ag"" = 0.812). 

The redshift distribution of the sources plays a vital 
role in modeling the observed angular correlations. Figure 
[T] shows the normalized redshift distributions, M(z), which 
were used in these calculations for both the HerMES and H- 
ATLAS data sets. The estimates for the redshift distribution 
for the HerM E S dat a sets are consistant with those used by 
ICoorav et al.l (|2010l ) in their analysis. Several estimates of 
the redshift distribution of objects in the H- ATLAS survey 



have been made. First, lAmblard et al] (|2010l ) used colour- 
colour diagrams to estimate the redshift distribution; how- 
ever this method includes only a subset of the sources used in 
the angular correlation analysis and is known to be slightly 
biased towards higher redshift objects. Second, iDve et all 
l|2010l ) used optical counterparts of sub-mm sources to esti- 
mate the redshift distribution; but again this method only 
uses a subset of all the objects and is biased towards lower 
redshift objects. So far, the best estimates of the redshift 
distribut ion, in that they u se all the available sources, are 
those bv lEales et al.l (|2010al ) who use Spectral Energy Dis- 
tribution (SED) fitting t o find a best fit reds hift distribution. 
We use the estimates bv lEales et all \2QIQ^ ) in this analysis, 
but the values for N{z) used here are for illustrative pur- 
poses only and should not be regarded as definitive given 
the preliminary state of the data. We have tried a variety 
of alternative models and find that, for reasonable choices, 
the results for tJobs are not especially sensitive to a partic- 
ular mean redshift in Af{z); this is probably because of the 
relatively slow evolution of the power spectrum at low red- 
shift in the concordance cosmology. On the other hand the 
results are sensitive to the width of the distribution in z; the 
wider the distribution over z the lower the amplitude of the 
angular correlation function. This is consistent with what 
you would expect if you consider that the clustering signal 
is more concentrated in narrower redshift bands. 

The Hersche l Multi-tiered Extragalactic Survey 
(|01iver et al.l |2010| . HerMES) will cover approximately 70 
deg'^ of the sky in three different wavelength bands. The 
Lockman-SWIRE field, which is one of the shallower fields, 
was covered by the initial science run and results of the 
cluste ring have recently been published in ICoorav et al.l 
(|2010t ) along with a Halo model analysis. Here we show just 
three of the available results, in Figure [2l to compare with 
the theoretical models we have discussed above: S250 > 
35mJy, S350/S250 > 0.85 and S350/S250 < 0.85. From the 
observations we see that the angular correlation has similar 
amplitudes in each of the three data sets; the middle panel 
has the highest amplitude and this most likely because the 
redshift distribution of the component sources is the most 
narrow of the three data sets making the clustering signal 
more concentrated and appear stronger. 

In Figure [2] four models are compared to the observa- 
tions: the short dashed line shows a simple power law fit, 
the dotted line shows the fit from the Halo Model, the solid 
line shows the non-linear approximation using the Fitting 
Function approach, and the long dashed line shows just the 
linear spectrum. In the top two plots, the power-law fits the 
data well at low angular scales, but under estimates at mid 
scales and over estimates at higher scales. This is similarly 
the case for the last plot but it also overestimates clustering 
on small scales. The Halo model provides a better fit; the 
low and intermediate angular scales both fit well, although 
in the top two data sets the model is slightly overestimating 
at large angular scales. The fitting function approximation 
also fits the data well. The fit at intermediate scales is per- 
haps not quite as good as the Halo Model, but then this 
approach has fewer free parameters; the differences in the 
behavior at large angular scales is a consequence of slightly 
different best-fit values for the bias parameter. However, an 
important point to note here is that the linear angular cor- 
relation function also provides a pretty good fit to each of 
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Table 1. Best fit parameters for the Halo and Fitting Function 
models for the HerMES results. Mmin is the minimum halo mass 
for which galaxies can form, Msat is the average mass of a halo 
with one satellite galaxy and as is the slope of the first moment 
of the satellite galaxy HOD. All masses are in units MQh~^. We 
see that the Mmin found for each data set are consistent between 
models, given an accuracy of 10^"'^ MQh~^ on the masses. 
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Figure 2. These plots show the angular correlation function 
'^obs (S) measured in the HerMES survey in the Lockman-SWIRE 
field. Sources were divided into three sets by flux density: S250 > 
35mJy (top), S350/S250 > 0.85 (middle) and S350/S250 < 0.85 
(bottom). Data is compared to four different theoretical models; 
the non-linear approximation using the Fitting Function (solid 
line), a power law (short dashed line), the Halo Model (dotted 
line), and the linear approximation (long dashed line). The Halo 
and Fitting Function methods fit the data well, as does the linear 
approximation. The power law over estimates at large scales and 
underestimates at mid angular scales. 



the data sets. By comparing the linear and non-linear mod- 
els we see that the transition between linear and non-linear 
regimes does not become evident until quite small scales (~ 
1 arcmin). The limited resolution of the Herschel telescope 
makes it difficult to probe the clustering regime on scales 
much smaller than this. 

In these models the amplitude of large-scale cluster- 
ing depends on the bias and that, in turns, boils down to 
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Figure 3. These plots show the angular correlation function 
iOohs W measured in the science demonstration phase of the H- 
ATLAS survey. Sources were divided into sets by flux density: 
S250 > 33m Jy (top) and S350 > 36m Jy (bottom). The data are 
compared to the Fitting Function model (solid line). 



the minimum halo mass through equation [9l We can there- 
fore use the measured clustering amplitude to get a rough 
estimate of the mass of halos hosting these galaxies. The 
values of the minimum halo mass corresponding to the the- 
oretical models in Figure [2] are shown in Table [1] We see 
that for each of the data sets, the minimum halo masses 
found from both the Fitting Function and Halo models are 
in agreement (within the errors) and suggest a value of 
Mmin ~ 10"-^=^° *MQft"^ The results using our Halo model 
ar e consistent with tho se found in the Halo Model analysis 
bv lCoorav etall (|2010l ). 
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The other results we have analyzed are from the Her- 
schel ATLAS survey l|Eales et al.ll2010l . H- ATLAS) which 
will cover 550 deg^ of the sky, in 5 different wavebands cov- 
ering the far infrared to sub-millimetre. The science demon- 
stration phase of the survey has just been completed (which 
covers roughly 1/30 of the final H-ATLAS data-set) and 
iMaddox et all (|2010l ') have released measurements of the 
angular correlation function of the galaxies observed. Fig- 
ure |3] shows different flux cuts for S250 > 33mJy (top) and 
S350 > 36m Jy (bottom). The results are quite noisy so none 
of the models fit particularly well - a detailed analysis will 
have to wait until the completion of the survey. This noise 
is thought to be due to galactic cirrus in the H-ATLAS data 
as the number of detections are very similar to the HerMES 
data sets and both surveys are similarly confusion limited. 
In this case, therefore we have just plotted the Fitting Func- 
tion results against the data. It was used to find approximate 
values of the corresponding minimum halo masses, which are 
10^^ '^ and 10^^-^ M^h'^ for the examples plotted in Figure 
[3]top to bottom respectively. The scale of the ojobs results for 
the S250 > 33mJy flux cut is significantly smaller that those 
in the other examples. This is predominately due to the 
wide bimodal redshift distribution (see Figure [1} although 
it does still suggest a slightly lower minimum mass ~ 10^^ ' 
Mq/i"^. This minimum mass is co nsistent wi t h tha t found 
in the cross correlation analysis bv lGuo et all (|2010ll . 



4 DISCUSSION AND CONCLUSIONS 

In this paper, we have compared models of the angular cor- 
relation function to data from the science demonstration 
phase of Herschel. We highlight a Fitting Function method 
which provides an improved fit to the data than a power 
law, and similar to that of the Halo model. It has just one 
free parameter, the minimum halo mass, compared to the 
two and three for a simple power-law and the halo model 
respectively. The halo mass is more meaningful physically 
as a parameter than those involved in the power law fit, 
so the Fitting Function is a much better method than the 
power-law for a quick-and-simple analysis. Although neither 
as sophisticated nor as flexible as the Halo model, it remains 
a useful tool that is perfectly adequate for modelling cur- 
rently available data. For example, the minimum halo mass 
found using our Fitting Function model is consistent with 
that found using the Halo Model. 

Another point of interest is that, in fact, the linear angu- 
lar correlation function also provides a reasonable fit to the 
data. Owing to the limited resolution of the Herschel tele- 
scope it is difficult to identify pairs of galaxies sufficiently 
close together on the sky to probe anything but the mildly 
non-linear regime. The currently available data provide some 
evidence of a transition between the linear and non-linear 
regimes, but they provide no unambiguous detection of the 
change of shape in a;obs that the Halo Model predicts. This 
does not mean the Halo Model is incorrect, of course, but 
what it does mean is that, at least for the time being, the 
paraphernalia involved in modelling intra-halo correlations 
is rather superfluous for these objects; simpler models can 
yield perfectly adequate results. 

Finally, we stress that the data sets to which we have ap- 
plied these models are preliminary. The biggest stumbling- 



block to a more complete analysis relates to the consider- 
able uncertainties in the redshift distribution N[z) of the 
sources involved. The choices we adopted for this analysis 
are for illustration only so the results should not be regarded 
as deflnitive. Further data, especially ancillary data provid- 
ing spectroscopic redshifts, will be needed before the precise 
nature of these galaxies can be determined. 
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